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O . Extending the analysis in our previous paper, we construct the entangle- 

5-H . ment thermodynamics for a massless scalar field on the Schwarzschild space- 

time. Contrary to the flat case, the entanglement energy E^nt turns out to 
be proportional to area radius of the boundary if it is near the horizon. This 
peculiar behavior of E^nt can be understood by the red-shift effect caused 
by the curved background. Combined with the behavior of the entanglement 
entropy, this result yields, quite surprisingly, the entanglement thermodynam- 
ics of the same structure as the black hole thermodynamics. On the basis of 
these results, we discuss the relevance of the concept of entanglement as the 

microscopic origin of the black hole thermodynamics. 
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I. INTRODUCTION 

One of the central problems in black hole physics is the identification of the microscopic 
origin of black hole entropy, which obeys the relation 



Sbh = j^A, (1.1) 

where A is the area of the event horizon and l-p\ := \G is the Planck length [|T|. (Hereafter 
we will set c = h = 1.) There are two important facts which suggest that the black hole 
entropy may have some microscopic origin. One is Hawking's argument [Q on quantum 
fields in a black hole spacetime showing that the black hole emits thermal radiation with 
the temperature Tbh which is proportional to that determined by the classical 1st law for 
black holes with a universal coefficient. In fact, the coefficient on the right-hand side of 
Eq.( |l.l|) is chosen so that Tbh coincides with the latter. For derivation of Tbh^ see [|,|1 as 
well as [0]. The other is the fact that the free energy calculated from the Euclidean path 
integral for the pure gravity system gives in the saddle point approximation exactly the 
same expressions for the temperature and the entropy as those given above. On the basis 
of these facts and their consistency, there have been proposed various candidates for the 
microscopic origin of the black hole entropy [p. Among them the simplest one is the idea 
of the entanglement entropy 0J^, which focuses on the entropy associated with states of 
quantum fields on black hole spacetimes and is strongly motivated by the first of the above 
two facts. 

The entanglement entropy itself is quite a general concept which is nothing but coarse 
graining entropy for a quantum system caused by an observer's partial ignorance of the 
information on the state. Now the idea is that the event horizon may play the role of the 
boundary of perception for an asymptotic observer so that its existence may give rise to the 
entropy. Indeed for simple models on a fiat spacetime [|^,|^ explicit estimations showed that 
the entanglement entropy Sent is always proportional to the area A of the boundary between 
two regions of a spatial section of spacetime: 



Sent = %A (1-2) 

where a is a cut-off scale for regularization and Cs is a model-dependent coefficient of order 
unity. Thus there is a clear similarity between Sent and Sbh- 

In order to see whether there is something deeper in this similarity, we have constructed a 
kind of thermodynamics of a space boundary, which we call entanglement thermodynamics, 
for a massless scalar field on fiat spacetime in our previous paper 0. There, by giving a 
suitable definition of the entanglement energy Eent, we determined the entanglement tem- 
perature Tent by imposing the first law of thermodynamics. Then we compared the system 
of thermodynamics quantities {Eent, Sent, T^nt) obtained by this procedure with the corre- 
sponding one for black hole, {Ebh = M, Sbh,Tbh)- This comparison showed that the 
entanglement thermodynamics on a fiat background possesses a totally different structure 
compared with the black hole one. 

It is not difficult to understood why such discrepancy occurs. Since we wanted to con- 
struct thermodynamics of a space boundary, we have defined the entanglement energy E^nt 
so that it depends only on quantum degree of freedom around the space boundary. As 
a result, Eent became proportional to the boundary area A unlike the black hole energy 
Ebh = M (X \fA. 

Since the dependence of the thermo dynamical quantities on the boundary area come from 
the very nature of the idea of entanglement, it may appear that the idea of entanglement 
cannot have any relevance to the black hole thermodynamics beyond the similarity of the 
expressions for entropy. However, it is not the case. It is because gravity is not taken into 
account in the above argument. In fact there is a good reason to expect that the inclusion 
of gravity improves the situation drastically. 

In black hole spacetime the energy of quantum field gets gravitational corrections, which 
depend on the definition of the energy. Taking into account of the fact that Eent depends 
only on modes around the boundary, the entanglement energy estimated on fiat background, 
Eent fiat, should bc identified with the energy Eent,h measured by a Killing observer located 



near the horizon. On the other hand, since the black hole mass used as the energy in the 
black hole thermodynamics corresponds to the energy measured by an observer at spatial 
infinity, it is natural to use the corresponding quantity Eent,oo- Due to the gravitational 
redshift, these energies are related by 

Eent,oo = {~9tt)h Eent,h — {~9tt)h -^ent,flat; (1-3) 

where {—gtt)\! is the well-known red-shift factor [||] for a signal emitted at the boundary. 

Let us assume that the boundary is at the proper distance a ~ /pi from the horizon. 
Then for the Schwarzschild metric we obtain {—gtt)h — ci/{'^M) oc l/yA. Hence from 
Eentfint OC A, if follows that Eent,oo OC \CA. Ou the other hand. Sent is independent of the 
position of an observer once a quantum state is fixed, so Sent oc ^ as before. Thus it is 
expected that the gravity effect modifies the structure of the entanglement thermodynamics 
so that it becomes identical to that of the black hole thermodynamics. 

On the basis of these observations, in this paper we construct the entanglement ther- 
modynamics for a massless scalar field on the Schwarzschild background, and compare its 
structure with that of the black hole thermodynamics. 

The paper is organized as follows. In the next section we describe the model used in the 
paper in detail, and clarify the definitions of the entanglement entropy and energy and the 
basic assumptions with brief explanations of their motivations. Then, after explaining our 
regularization scheme, we derive general formulas for the entanglement entropy and energy 
of the regularized system and estimate them numerically in §3 and §4. On the basis of these 
results, we compare the structure of the entanglement thermodynamics for the massless 
scalar field in the Schwarzschild background with that of the black hole in §5. 

II. MODEL CONSTRUCTION 

In order to construct the entanglement thermodynamics of a black hole, we must make 
clear how to implement the idea of entanglement into a quantum field system on a black 



hole spacetime, and give the definitions of the basic thermodynamical quantities in terms of 
the quantum fields. 

A. Model description 

The basic idea of the entanglement thermodynamics is to express the thermodynamical 
quantities for a black hole in terms of expectation values of quantum operators dependent 
on the spacetime division as in the statistical mechanics modeling of the thermodynamics 
for ordinary systems. Therefore we must specify how to divide spacetime into two regions 
and with respect to what kind of state the expectation values are taken. 

According to the original idea of entanglement, it is clearly most natural to consider a 
dynamical spacetime describing black hole formation from a nearly flat spacetime in the past 
infinity, and divide the spacetime into the regions inside and outside the horizon. In this 
situation, if we start from the asymptotic Minkowski vacuum in the past, the entanglement 
entropy associated with the division of spacetime by the horizon acquires a clear physical 
meaning. However, this ideal modeling cause difficulties. The most serious one is caused 
by the occurrence of the Hawking radiation: its contribution to the entanglement entropy 
diverges. Of course, if the backreaction effect is properly taken into account, this divergence 
may disappear. But such a modeling is intractable at present and is beyond the scope of 
the present paper. 

In order to avoid this difficulty, we replace the problem by a stationary one on the basis 
of the following arguments. First, consider a situation in which a thin spherical shell with 
a tiny mass 5M infalls toward the origin to form a Schwarzschild black hole with mass 5M. 
In this situation, the timelike Killing surface S with r = 2GAM in the Minkowski region 
turns to the horizon after the passage of the shell. Two things happen in this process: a 
division of spacetime by the horizon and a change of geometry. Among these two things, 
the latter is responsible for the Hawking radition, because the Hawking radiation is essen- 
tially determined by the geometry near the horizon in the geometrical optics approximation. 



Meanwhile, the division of spacetime produces the entropy of the black hole in our approach 
based on the entanglement. Thus, we can determine the entropy of the black hole avoiding 
the above difficulty caused by the Hawking radition, if we calculate the entanglement en- 
tropy associated with the division of Minkowski spacetime by the timelike surface S starting 
from the Minkowski vacuum. Because we need to introduce a cut-off length a of the order 
of the Planck length /p;, a natural minimum value of AM is determined by the condition 
r^ ~ a^. The corresponding value of the entropy becomes of order unity. 

Next, let us consider a thin spherical shell with a tiny mass AM infalling toward the 
center in a Schwarzschild background spacetime with mass M. After the passege of the 
mass shell, a Schwarzschild black hole with mass M -|- AM is formed. As in the previous 
case, this process produces a new horizon which divides the spacetime by a null surface 
corresponding to the stationary timelike surface S with r = 2G{M + 5M) before the passege 
of the shell, and slightly modiffes the geometry. If we apply the reasoning in the previous 
case to the present situation, we can expect that the entropy of the new black hole is given 
by the entanglement entropy for quantum ffelds in the Killing vacuum state, associated with 
the division of the Schwarzschild spacetime by the timelike surface S before the passage of 
the shell. Here, the Killing vacuum state should be taken in order to avoid the Hawking 
radiation. Now the minimum value of AM is determined by the condition that the increase 
of the area of the horizon is of order a^, which corresponds to an increase of entropy of order 
unity. For this minimum value, the proper distance of S to the original horizon becomes of 
order a if a ~ /p;. 

On the basis of these arguments, in this paper, we consider a quantum ffeld in the Killing 
vacuum state in the Schwarzschild background with mass M, 

d-s^ = -(1 - ro/r)dt^ + (1 - ro/r)-^dr^ + r^dn^, (2.1) 

where tq = 2GM, and calculate the entanglement entropy and the entanglement energy 
to be deffned later associated with the division of the spacetime by a timelike surface S 
at a proper distance of the order of a cut-off length(~ Ipi) to the horizon. Since there is 
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no definite criterion on the exact position of tlie boundary in our framework, we will also 
investigate the influence of the variation of its position. 

Here a subtlety occurs: if we require that the quantum field is in the Killing vacuum state 
on the whole extended Kruscal spacetime and take the horizon as the boundary surface, then 
the entanglement entropy vanishes because the state is expressed as the tensor product of 
the Killing states in the regions I and II ( Figure^). This triviality comes from the fact that 
the quantum degree of freedom on the horizon is completely eliminated for such a state. In 
order to avoid this, we restrict the spacetime into the region I, and replace the boundary 
by a time-like surface S at a proper distance of the order of the cut-off length (typically, 
Planck length) of the theory to the horizon [} The boundary S can be regarded as a would-be 
horizon when matter fields (a real scalar field in our analysis in this paper), which contribute 
to the entanglement entropy, falls into the black hole. Evidently, the would-be horizon must 
be very close to the horizon since our analysis is semi-classical one. However, because of the 
existence of the cut-off length there is the minimum proper length between the would-be 
and the present horizon: it must be of order of the cut-off length. Hence it is natural to take 
as the boundary a timelike surface at a proper distance of the order of the cut-off length to 
the horizon. This prescription is expected to give a correct estimate of the thermodynamical 
quantities for the case in which the spacetime is divided by the horizon and the quantum 
state inside the horizon is arbitrary specified because the entanglement quantities depend 
only on the degrees of freedom within the distance of the order of the cut-off length from 
the boundary. However, we should keep in mind that we have no definite criterion regarding 
the exact position of the boundary. To minimize this ambiguity, we will also investigate the 
influence of the variation of the boundary position. 



^ Our present model clears such a criticism |jl^ for works done before ||7|-p|,p!l|| that in these works 
the boundary had always been chosen as a timelike surface, contrary to the event horizon, which 
is a null surface. 



As a matter content we consider a real scalar field described by 

S = -lj{d'<l^d,<jy + fi'},/^dx' . (2.2) 

The mass fi does not play an essential role since a typical length scale controlling the 
entanglement thermodynamics is much smaller than the Compton length of a usual field. 
Therefore we just set /i = in the numerical computations. 

B. Entanglement entropy and energy 

There is no difficulty in generalizing the entanglement entropy for the flat case ^^ 
to the case of Schwarzshild background. In contrast, with regards to the deflnition of the 
entanglement energy, we cannot simply extend the deflnition used for the flat case [Q due 
to the difference in the method of implementation of the idea of entanglement as well as due 
to the existence of the black hole mass in the present case. 

Before discussing the deflnition of the entanglement energy, let us flrst recall the standard 
general procedure of building a statistical mechanics model for thermodynamics. Let A be 
a subsystem surrounded by the environment B. For a quasi-static process, the flrst law of 
thermodynamics for the subsystem A is expressed as 

dE = TdS + fdx . (2.3) 

Here E, S and T are the thermodynamical energy, entropy and temperature of A, respec- 
tively, / is a generalized force applied to A, and x is a generalized displacement. Now the 
statistical mechanics model of this thermodynamics is constructed by the following sequence 
of steps: 

(a) Assign a microscopic model to A and construct a Hamiltonian system for it. 

(b) Make statistical assumptions (e.g. the principle of equi-weight) which deflne averaging 

procedures for the microscopic variables. 



(c) Identify the thermodynamical variables with quantities constructed from the micro- 

scopic variables introduced in (a) along with the averaging procedure in (6). 

(d) Reproduce the relations between the thermodynamical variables (thermodynamical laws 

in particular). 

Since the basic idea of the entanglement thermodynamics is to construct a statistical 
mechanics model of the black hole thermodynamics using the idea of entanglement, we 
should also follow a similar procedure to this general one in defining Ef,nt- 

In our case it is natural to regard the scalar field in the region inside the boundary E 
as the subsystem A, and the scalar field in the region outside S as the environment B. 
Then, if we decompose the Hamiltonian H corresponding to the Killing energy to the parts 
dependent on the microscopic degrees of freedom inside S, those outside S and both as 
H = Hin + Hout + Hinti the most natural choice of the microscopic variable for Eent is the 
operator i7j„. Clearly the averaging procedure is given by taking the expectation value with 
respect to the Killing vacuum of the total Hamiltonian H. Hence we are led to the definition 

where : Hin '■ denotes the normal ordering with respect to the ground state of Hin- 
Here note that the choice 

(F) E^JJ = (: H^^t :) 

is essentially same as (/) in spite of its appearance, as we will see in Sec.|^. It is because 
Eent measures a sort of disturbance to the vacuum state caused by the boundary S, and 
its value is determined just by the modes around S. Hence the difference between (J) and 
(/') comes from a tiny difference in the redshift effect on the modes just inside S and just 
outside S. 

These choices should be contrasted with the definition of the energy in the standard 
argument of black hole thermodynamics. There the energy is given by the black hole mass 



M, which is quite natural in the framework in which the backreaction of matter on gravity 
is consistently taken into account. However, in the present case, the option 

(II) E^IP = M, 

lacks charm in our framework because it is not related to microscopic degrees of freedom. 
Of course, it may be reasonable to replace the definition (I) by 

(III) E^^f =M+(:i7,„:), 

in the sense that it represents a kind of total energy of the system consisting of the scalar 
field and the gravitational field. 

Here we should comment on a subtlety of the role of M in our model. Since we are 
neglecting the backreaction of the matter on gravity, the black hole mass M determining 
the background geometry gab might be regarded as the external parameter x in the above 
expression for the 1st law of thermodynamics. However, since the Killing vacuum state 
depends on M as we will see later, M becomes a function of thermodynamical variables Sent 
and E(.nt- Thus in effect the term fdx in Eg. ( ^.31 ) does not appear in the present case ||12[|.p| 



Because there is no firm ground to pick up one of these, we will investigate all of them. 
Further, for the sake of comparison with the result in the flat spacetime models discussed 
in our previous paper |P], we also calculate the quantities defined by 

(IVi) e'^Z'^ := (: H :),,, 

(IV2) E^Z'^ = {:H,n:) + {:H,ut:). 

Here : H : and : Hout denote the normal products with respect to the ground states of H 
and Hout^ respectively, and {■) pi indicates the expectation value with respect to the density 



^ It is not a surprise that an external parameter becomes a function of variables describing 
thermodynamics state. For example, in the Rindler spacetime, the acceleration a is determined by 
the Rindler temperature. 
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matrix p' := pin® Pout, where pj„ and pout are the reduced density matrices for the inside 
and the outside of S, respectively 0. 

Finally we comment on the regularization. In dealing with the matter field, we naturally 
encounter the ultraviolet divergence. We adopt here the cut-off regularization by introducing 
the length scale a, which is supposed to be of order of l-p\. On dimensional grounds we 
expect that Sent oc A/a^ and Eent oc ^/A/a'^. Thus the entanglement temperature T^nt = 
dEent/dSent IS cxpccted to be independent of the cut-off scale a. We will see later that this 
is the case except for the definitions (//) and (///)• 

III. ENTANGLEMENT ENTROPY 

A. Basic formulas 

In order to estimate the entanglement entropy and the entanglement energy for the scalar 
field on the Schwarzschild spacetime, we first regularize the field theory described by Eg. ([2. 2]) 
and reduce it to a discrete canonical system described by a Hamiltonian of the following 
form : 

N . -, 

^0= E ^^^''PAPB + -^ABgV, (3.1) 

A,B=1 ^^ ^ 

where {{q'^,PA)} {A = 1,2,- ■ ■ ,N) are canonical pairs. Here both g'^'s and p^'s are of 
physical dimension [L°], and the parameter a is the cut-off length of order of lp\. 

Because of the spherical symmetry of the system, if we expand the scalar field in terms 
of the real spherical harmonics as 

0(t, p, e,ip) = Y^ <Pim{t, p)Zi^{e, if), (3.2) 

the action becomes a simple sum of the contributions from each mode (f)im- Here Zim = 
^/2^Yim, V^'^Yim for m > and m < 0, respectively, with Ziq = Yiq, and p is a suitable 
radial coordinate. Hence if we discretize this radial coordinate, we obtain a regularized 

system. 
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As the radial coordinate, we adopt the proper length from the horizon, which is related 
to the circumferential radius coordinate r by 



where v := J\ — ■^. In this coordinate the metric Eq.( |2.1| ) is written as 



2 



ds' = -v''dt^^dp'^rl- — , (3.4) 



'\-v 



where v is understood as a function of p through Eq. (|3 
Plugging Eqs.(U) and Q into Eq.Q, we get 



S = -J2 j dt dp 

Im 



rUL i rlv ,^ , x2 , (ui , IX , ('^o/i) 



^(1 _ ^ij^ - \(l_^2)2 ^-P^^rn, ^ y,^ ^ x; ^ (1-^2)2^ ^^L 



{d,<Pimr +/(/ + !) + yf^^^ ^0: 



(3.5) 



Comparing this equation with Eq. (|3.1|) we see that the most suitable configuration variable 
is the dimensionless one defined by 

Hereafter we adopt a boundary condition that ipim{t,p) is finite at the horizon, which corre- 
sponds to the Dirichlet boundary condition for the original variable (f)im{t, p)- To be precise, 
in defining Hamiltonian for our system, which must be self-adjoint, there are two options 
for the boundary condition of (pimit, p) at the horizon: the Dirichlet boundary condition and 
the Neumann boundary condition. In our case we adopt the former since the later makes 
the Killing energy divergent. 

Now it is straightforward to perform the canonical transformation. The Poisson bracket 
relations become 

bPlm{p), T^l'm'ip')} = Sll'6mm'S{p - p) , 

{^Zm(p),^i'm'(p')} = ) 
{T^lm{p),'^Um'{p')} =^ , (3.7) 
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where Tiim is a momentum conjugate to ipim with dimension [L ^]. In terms of these canonical 
quantities the Hamiltonian is expressed as 



H = j:Hi 



Im 



Im 



Hlm= 2 y dpTrfrnip) + ^ / dpdp'^lm{p)Vlm{p, p')^lm{p) 



(3. 



where 



i'lm{p)Vlm{p, P')i^lrn{p') = S{p, p') 

,2 



:{d,{u'/\l-u')^,Jiy 



1 - Z/2)2 

+ ^[l(l + l)(l-u')' + {ropY}t 



(3.9) 



We regularize this system by replacing it by a difference system with respect to p with 
spacing a. To be precise, we make the following replacements: 

p -^ Aa, 



tPiUp =iA- l/2)a) - 


^qL, 


Tllra{p ={A- l/2)a) - 


-^ PlmA, 


u{p={A-l/2)a)- 


^ l^A, 


6{p = Aa, p = Ba) - 


^ Sab /a, 



where A runs over the positive integers. To achieve a better precision, we adopt the middle- 
point prescription in discretizing the terms including a derivative: we replace a term, say, 
f{p)dpg{p) by /a+i/2 ■ ^{dA+i — 9 a)- Further, in order to make the degrees of the system 
finite, we impose the boundary conditions q^^^^ = 0. In the numerical calculation A^ is taken 
to be sufficiently large so that this artificial boundary condition, which is required just for 
a technical reason, does not affect the results. 

In this manner we get the hamiltonian in the desired form: 



Im 

N 

A,B=1 



(Im) 



^0 PlmAPlmB + ^ ^AB QlmQlm 



(3.10) 
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where 



N 



N 






A,B=1 



A=l 



^A+1/2 



^^ ^A+1/2) 



1/2 



1/2. 



,2\^A 



^A;^(i-^i-M)g."j"^-^r(i-^i)g^ 



+ (^j z.l(/(/ + l)(l-z.l)^ + (ro/i)>/ 



2A.A2 
m 



Here \/('™) becomes the positive definite, symmetric matrix whose exphcit form is 



Ai S(') A2 



2a 



(l/jM) = f!^ 



\ 



Aa-1 S^ Aa 



.(0 



2a2 



z/a(1 - z^: 



2\2 



^A+1/2 



+ 



^^^-1/2 



;i - ^1+1/2)^ (1 - ^1-1/2)^ 



+ iz.l(/(/ + l)(l-z.l)2 + (ro/.f) 



A 



A 



2a2 



:^'^A+l/2(«^A+l«^AJ 7^ -2 7^ • 



(1 - 4+1/2)^ 



(3.11) 



Here note that the mass term in S^ is neghgible compared with the first term iiafi « 1. 
Therefore we will simply set /i = in the numerical calculations. 

B. Formulas for the entanglement entropy 



In the discretized system^ {{q^iPa)] Hq\ (A = 1, 2, ■ ■ ■ , A^) with Eq. (|3.1CI|) , the position 
of the space boundary S is given by p = p^ := ^sa with ns = 0(1), which the set of modes 
{q^} {A = 1,- ■ ■ ,N) into two subsets, {g"} (a = 1, ■ ■ ■ , ub) and {g"} {a = ub + I,- ■ ■ , N). 
Here we regard {q°'} and {g"} as the inside modes and the outside modes, respectively. 

The Killing vacuum for the continuous system corresponds to the ground state of the 
Hamiltonian Hq in this system. Hence its density matrix is given by 



For the notational simplicity we will often omit the suffices {l,m) if no confusion occurs. 
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iB 



P(K},{r})= det- exp 



iy\V2 



vr 



-^w^AB(gV + g'V'') 



(3.12) 



where W = {aVY^'^. In accordance with the sphtting of {q^} into {q"^} and {(?"}, the 
matrices F, W and its inverse PF"^ naturally split into four blocks as 



{V) 



AB 



V 



K 



(1) 



ab 



i^int)al 



\ 



{V^nt') 



ab 



K 



(2) 



■P ) 



{W) 



{w-^ 



AB 



AB 



Aah B, 



ab ^af3 



\ 



\ {B )ab Dap I 
( 



(3.13) 



V 






\ 



-.A\ f^fB 



Taking the partial trace of p({g }, {q' }) for the inside modes {q"^} (a = 1, ■ ■ ■ , ra^), we get 



the reduced density matrix 



^\-l/2 



Pred({g"},{g'^}) = (det7r^)"' exp 



-liD-'Uiq-q^ + q'W^ 



a I ii^jP^ 



--{B^A-'BUiq-qTiq-qr 



Now the entanglement entropy associated with the boundary S, Sent '■= — Trpredlnpred, 
is given as follows |0,|[. Let {A^} {i = 1,- ■ ■ ,N — ub) be the eigenvalues of a positive definite 
symmetric matrixFI A, 



A:=D'/^B^A-'BD'/\ 



(3.14) 



Then it is easily shown that modes labeled by (/, m) contribute to the entangle entropy by 
the amount 



^ The corresponding expression in ref. ("A'^t := {M^^)''''Ncb") reads A = DB'^A^^B in the 
present notation. This definition does not give a symmetric matrix and should be replaced by 

"A«'' := (M-i/2)'^^Ar^rf(M-V2)df"'^ namely Eq.(|3l|). 
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N-riB 

'11 



c(0 _ Y^ c 



Si = -- — ^In/Xj -ln(l -/Xj), (3.15) 

1 -/ii 

where \ii := \^^ {\J\ + Aj — Ij . (Note that < /ij < 1.) Clearly 5'eJt is independent of 
m = {—I, — / + 1, ■ ■ ■ , / — 1, /). The entanglement entropy is given by 

oo 
Sent = T.i^l+l)Sii. (3.16) 

/=0 

From Eq. (|3.11|) , one can easily show that 

Sent ^ O ({la/ro)~'^\n{la/ro)) ^^ /a/ro ^ oo. 



Thus the infinite series Eg. ( p. 151 ) actually converge so that we can safely truncate them at 



some appropriate /, depending on the accuracy we require and the ratio ro/a we set. 

C. Numerical estimation 

Using these formulas, we have evaluated Sent numerically and have examined its de- 
pendence on the area of the boundary, A = Anrj^. In this calculation the outer numerical 
boundary is set at A^ = 60. The summation in / in Eg. ( p. 161) is taken has up to / = [lOro/a] 



([ ] is the Gauss symbol). From the above asymptotic behavior of S^.^^, this guarantees the 
accuracy of 10%. 

The result is shown in Figure ^j. From this figure we see that Sent is proportional to 
A/a"^ if we change tq with fixed ub, and its coefficient only has a weak dependence on ub- 
In particular, for the limit n^ = 1, we get 

Sent ^ 0.024^/0^ . (3.17) 

This result is essentially the same as our previous result for models in a fiat spacetime 
including the numerical coefficient 0. This can be understood in the following way. 
Let us make a coordinate change from r to x defined by 
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r 
To 



(x + 1) 
4a; 



or 



2r //2r \^ 

a;= l + W II -1. 

ro V Vro 



Then Eg . (12.11 ) is rewritten as 



ds^ 



Jh -L \ 1 , o o / -L I t// 



x+ 1 



rft' + r^ 



2x 



{dx' + x'dn^) 



Note that r = 0, ro and oo correspond to a; = 0, 1 and oo, respectively. It is easy to see that 
the Hamiltonian (Eq.( p.8|) ) is given in this coordinate system as 



H = 22 Him, 

Im 



H, 



Im 



d^ 



6Ax*{x - \] 



'.pf 



Im 



1 ( {x+lf ' 

2 \ Ax 



{d^ip 



Im) 



/(/ + !) 

e 



-Vlr. 



(3.18) 



— 4r"tin 0/m and ipim ■■= ro0/m in 



where ^ := rox, and Pim and v?im are expressed as Pim ■- 64x4(a;-i) 
terms of (pim in Eq.( p.2|) . Here note that the vacuum state is only weakly dependent on the 
prefactor Z+tV' ^^ Eq. (|3.18|) . If we neglect this prefactor, the vacuum state is determined 
by the Hamiltonian which coincides with that for the flat spacetime at x = 1. On the other 
hand, Sent depends on the modes in a thin layer around the boundary S, whose typical 
thickness is a few times of a ~ /pi). Therefore, when S is near the horizon, the value of Sent 
should be well approximated the fiat spacetime value. 



IV. ENTANGLEMENT ENERGY 



In this section we give formulas for the various definitions of the entanglement energy 



introduced in Sec. |IIB| , and estimate their values numerically. 

First we derive formulas for the entanglement energies corresponding to the definitions 

(I) and (r): 
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E. 



\- -tiin •/; 



(I) .. 
ent • 

Eent '■— V tiout '■) ■ 



(4.1) 
(4.2) 



By rescaling the variables \(l^\ in § pT 



as 



g^:=5^^(W^^/^)^^g^, 



the expression of the density matrix for the vacuum state Eq.( |3.12| ) gets simphfied as 



N 

n 

C=l 



r 1 
g-^} \P\ {?-'^}) = n --^/'exp \-\ [{ff + {q'^Y} 



and the normal ordered Haniiltonian : Hin : is represented as 



• J^ in • 



2a \dq 



«'S'A(i,+^"--'^v 



dqb 



^ ^%f 



dq 



% + -a?" 



Here w^^^ is the positive square-root of aV^^\ and the matrices U and w^^^ are defined as 

V /Ca V JbD ' 



D_B- 



Hence the matrix elements of : Hm : p with respect to the basis \q^) are expressed as 
[q^} I : H,^ : p| {g^}) = j- { [{w^^^ + l)Uiw^'^ - 1)] ^^ g^"^ + Tr [f/(l - ^(^))] } 



N 



X [| TT^^/^exp -(g 



^C\2 



From this we obtain^ 



Eil^t = l(^Udfy{q^}\:H.^,r.p\{f 



Aa 



[aV;;l\Ar + A^,5'^'-2w^^5'^\ 



(4.3) 



See Eq.(|3.13|) for the definitions of the matrices A, A, D and D. 



Similarly i?g„j is expressed as 



<n 



4a 



-'ent 



,(2)^^/3' 



where w^"^^ is the positive square-root of aV^'^\ 

Eent corresponding to the definition (///) are simply related to E^.^^ by 



(4.4) 



E. 



(HI) 



ent 



M + Ei'J, 



(4.5) 



Further, Eent corresponding to the definitions (IV) have already been given for the fiat case 
H]. They are expressed as 



E. 



{IVi) ._ 



ent 



E. 



im) ._ 



ent 



Tr [: Htot : p'] , 

Tr [(: Hi^: + : H,,t :)p] 



(4.6) 



Here Htot ■= Hi„ + Hout + Hint and p' := p^ 
they become [H] 



Pout 



In particular for a vacuum state, 



^ent 



E. 



ent 



It, 

2 
2 



VintB 



VLB 



— Tr 
2a 



r«)i/2l+Tr[(aV(2))i/2l -Tt[W] 



(4.7) 



where Vint, V^^''^^ and W are given in Eq.( p.l3| ). 

Like the entanglement entropy (Eq.( |3.16| )), the entanglement energy is also given by the 
summation of each contribution specified by (/,m): 



Eent = f:{2l + l)Eil. 

1=0 



(4.8) 



(11) 

ent 



Here E^nt represents any kind of the entanglement energy defined above (except for E, 

M and E^nt = M + E^nt) and E^nt is its (/,m)-contribution. From Eq. (|3.11| ), it is easily 

shown that 



^Note that E^nt^ and E),nt^ here are the generaUzations to the curved background case of E^nt 
and Elnt, respectively, in Ref. [^. In the same way the subscripts Hn'' and 'out' correspond to '1' 
and '2' in Ref. il. 
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^^Int ~ O [ila/ro) ^) as la/r^ -^ oo. 
Thus the summation with respect to / in (|4.8| ) converges. 



With the helps of these formulas, we have numerically evaluated -EeJ, -Ee„i E^^nt ^'^^ 
Eent^ ■ Now we have taken the numerical outer boundary at A^ = 100. The truncation in 
the /-summation is the same as for Sent (up to / = [lOro/a]), which implies that the accuracy 



is about 10% from the above asymptotic estimate for E)^, 



(0 



nt- 



The results of numerical calculations are shown in Figure ^, Figure ^, Figure ^ and Figure 
p. In these figures, rBE^nt is plotted as a function of {rs/ay for ub = 1, 2, 5. All of these 
figures show that rBEgnt is proportional to riBirB/aY'- 

4nl ~ 0.05(^5 - ll2)rBla\ 

Ei'J ~ O.ObiriB + l/2)rB/a^ 
e'^Z'^ ~ Q.2nBrBla\ 
E'^Z'^ ~ O.lnsr^/a^. 



(4.9) 



Note that the identity E^nt = -^Lt + E^J holds from the definitions. 



From these equations we immediately see that the values of E^ni and E^nt coincide 



except for a tiny difference independent of n^. As was mentioned in Sec. [ll B| , this difference 

is understood by the gravitational red-shift: E^Jt comes from the modes just inside E while 

(I') 
E^nt originates from the modes just outside S. In the present numerical calculations, it 

means that E^nt and E^nt are determined by the modes at p = {ub — l/2)a and p = (ub + 

l/2)a, respectively(see Eq. (|3.11|) ). Hence, taking account of the fact that the contribution of 



each mode to the entanglement energy is proportional to the red-shift factor at its location, 
the ratio of E^nt and E^nt should be approximately given by 

Ei'J, : EiQ ~ u{p = {ub - l/2)a) : u{p = {tib + l/2)a) 

~ (riB - 1/2) : (ub + 1/2). (4.10) 
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This is consistent with the above numerical result. 

This argument is also supported by the numerical result for the fiat spacetime model 
shown in Figure |^. In this figure the values of El,n{ and E^..^^ for a massless scalar field in 
the Minkowski spacetime with S = S*^ are plotted. In this case there is no gravitational 
red-shift effect, so we expect that E'en! = -^Lt ) ^^ confirmed by the numerical calculation. 

V. COMPARISON BETWEEN ENTANGLEMENT THERMODYNAMICS AND 

BLACK-HOLE THERMODYNAMICS 

Now on the basis of our results for Sent and Eent, let us compare the structure of the 
entanglement thermodynamics with that of black hole. 

From the numerical results in §3 and §4, the entanglement entropy Sent and the entan- 
glement energy Eent for our model are expressed as 

Sent = Cs(^-f] J- (C5^ 0.096), (5.1) 



tpl 



2 



rs ^ frB\ (hi\ ,.2 



E.„. ^ C.^ ^ C. ^-^J ^r^J Mc^ (5.2) 

where 

d^^ ~ 0.05, dp ~ 0.15, 

4'^^^ -0.4, cP'^c^O.2. (5.3) 

Here and hereafter we will only consider the horizon limit, ub = 1- 

It is helpful to keep in mind that the case Cs = Ce = 1 along with a = Ipi corresponds 
to the black hole thermodynamics.]^ 



^ Strictly speaking A in Eg. ( [5. 2D is the area of the boundary and it differs from tlie area of the 
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From these expressions and the first law of thermodynamics 

dE = TdS (5.4) 

the entanglement temperature Tent is determined as 

Ce 

Tent = -p^TsH- (5.5) 



C, 



Thus we get 



s 



T^il ~ 0.52TBH, T^nt^ ^ I.GTbh 



n'J'> ^ A.2Tbh, nut" ^ 2.1TBH- (5.6) 

These results have several interesting features. First of all we immediately see that 
the entanglement thermodynamics on the Schwarzschild spacetime show exactly the same 
behavior as the black hole thermodynamics, as summarized in ra6/e|l[ This behavior is just 
what we expected from the intuitive argument in the introduction: the gravitational redshift 
effect modifies the area dependence of Eent so as to make the entanglement thermodynamics 
behave just like the black hole thermodynamics. 

Second it should be noted that the temperature Tent becomes independent of the cut- 
off scale a only for the options type (J) and {IV) (i.e. (/), (/'), (IVi) and (/V2)). This 
indicates that (II) and (III) are not good definitions of the energy in our approach. It is 
also suggestive that the average |(TgJ + T^nt) gives almost the same value as Tbh- This 
averaging corresponds to averaging out the difference in the red-shift factors for the one- 
mesh 'inside' and the one-mesh 'outside' of the boundary. Therefore such an averaging may 
have some meaning. 



horizon. However we here ignore the difference since it is totally negligible as compared with A 
itself. Similarly ^'^~^'° = O ((7^)^) so that we can safely set r^ ~ tq in the present context. 
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To summarize, our model analysis strongly suggests a tight connection between the 
entanglement thermodynamics and the black hole thermodynamics. It is worth emphasizing 
that matter fields on a black hole spacetime capture the background characteristics. Of 
course, our model is too simple to give any definite conclusion based on it. In particular, the 
ambiguity in the definition of the energy comes from neglecting backreaction of the quantum 
field on gravity. Further, even in the fixed background framework, our model is too simple 
in that its thermodynamics is essentially controlled by one parameter corresponding to the 
black hole mass. It is obviously useful to see whether the entanglement thermodynamics is 
consistent with the black hole thermodynamics for models with more parameters, such as 
those on the Reisner-Nordstrom spacetime or on the Kerr spacetime before attacking the 
difficult task of going beyond the semi-classical approximation. 
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FIGURES 
FIG. 1. The Kruskal extension of the Schwarzschild spacetime. We consider only the region / 

(the shaded region). As the boundary S we take the hypersurface r = rs- 

FIG. 2. The numerical evaluations for Sent for the discretized theory of the scalar field in 
Schwarzschild spacetime. Sent for ub = 1,2,5 is shown as functions of A/a^, where A is area of 
the boundary. We have taken A^ = 100 and performed the summation over / up to lOro/a. 

FIG. 3. The numerical evaluations for -Egnt ^^ the discretized theory of the scalar field 
in Schwarzschild spacetime. vbE^.^^ for ub = 1,2,5 is shown as functions of {tb/cl)'^, where 
tb = r{p = nBd). We have taken N = 100 and performed the summation over / up to lOro/a. 

FIG. 4. The numerical evaluations for i?g„/ for the discretized theory of the scalar field 
in Schwarzschild spacetime. tbE^^I for n^ = 1,2,5 is shown as functions of {rB/o)"^-, where 
vb ^ r{p = nBd). We have taken N = 100 and performed the summation over / up to lOro/a. 

FIG. 5. The numerical evaluations for E),^^ for the discretized theory of the scalar field 
in Schwarzschild spacetime. tbE^.^^ for ub = 1,2,5 is shown as functions of (rB/a)"^, where 
rB = r{p = nBd). We have taken N = 100 and performed the summation over I up to lOrg/a. 

FIG. 6. The numerical evaluations for E),^i . rBE).^^ for ub = 1,2,5 is shown as functions 
of (rB/a)'^, where vb = r(p = ubo). We have taken N = 100 and performed the summation over 
I up to lOro/a. 

FIG. 7. The numerical evaluations for E^^J^ and -E'en/ i'^ ^^e Minkowski spacetime. aE^^^ and 
aE^^l are shown as functions of (rB/a)'^, where tb is the radius of the sphere which provides the 
division of region. We have taken A^ = 60 and performed the summation over I up to lOro/a. 
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TABLES 



TABLE I. Comparison of two kinds of thermodynamics 



Entanglement in Schwarzschild spacetime 



Black-Hole 



Varied 

Fixed 

S 

E 

T 



A 
a 

oc A 

oc ^-1/2 



A 

hi 

oc A 

OC^l/2 
oc ^-1/2 



Varied 

Fixed 

S 

E 

T 



a 
A 

oc a 
oc a~ 
oc a' 



-2 







Zpi 
A 



oc / 



PI 



oc / 



PI 



oc/Oj 



Varied 

Fixed 

S 

E 

T 



a 

Eent 

OC a^ 

oc a^ 

oc a^^ 



lp\ 
M 



oc /pi 



oc/^i 



oc / 



-2 
PI 
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